KfWCI  I#€  |oANI  (|  pUMM  RCSfMCH  INST 

^AVcRSilISu  SU^Jt  TO  8-?iV?Tu) 

,  .  SEP  82  R  C  OAVIDSOW,  K  T  TSAN8  N0001H*79>C-OS9S 

tfffcLASSMXfO  SM-254-82-2R3-LJ  N0°°l  79  C-0555 


SCIENCE  APPLICATIONS ,  INC. 


Plasma  Research  Institute 
934  Pearl  Street,  Boulder,  Colorado  80302 


82  09  94  062 


SAI-254-82-243-LJ 

PRI-49 


September  1982 


NONLINEAR  BERNSTEIN -GREENE -KRUSKAL  WAVE  EQUILIBRIA 
SUBJECT  TO  GLOBAL  ENERGY 
AND  MOMENTUM  CONSERVATION  CONSTRAINTS 


Ronald  C.  Davidson*  and  Kang  T.  Tsang 


Science  Applications,  Inc. 
Plasma  Research  Institute 
934  Pearl  Street 
Boulder,  Colorado  80302 


L 


*****’■♦«*»  PflVqftS'!  ‘  ' 


•Permanent  address: 


MIT  Plasma  Fusion  Center 
Cambridge,  Massachusetts  02139 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWiM  D.i.  Ent.r.d; 


j  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

SAI-254-82-243-LJ/PRI-49 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtltla) 

NONLINEAR  BERNSTEIN-GREENE-KRUSKAL  WAVE  EQUILI¬ 
BRIA  SUBJECT  TO  GLOBAL  ENERGY  AND  MOMENTUM 
CONSERVATION  CONSTRAINTS 

5.  TYPE  OF  REPORT  b  PCRIOO  COVEREO 

Technical  Report 

e.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfaJ 

Ronald  C.  Davidson  and  Kang  T.  Tsang 

s.  contract  or  grant  number;.; 

N00014-79-C-0555 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Science  Applications,  Inc. 

Plasma  Research  Institute 
_  934  Pearl  Street,  Boulder,  CO  80302 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  b  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

ONR  Physics  Program  Office 

800  N.  Quincy  Street 

Arlinqton .  VA  22217 

12.  REPORT  DATE 

September  1982 

13.  NUMBER  OF  PACES 

35 

14.  MONITORING  AGENCY  NAME  S  ADDRESS;//  dlUatant  /ram  Controlling  Oltlco) 

IS.  SECURITY  CLASS,  (ol  (hi.  roport) 

unclassified 

is*,  dcclassification/oowngrading 
schedule 

IS.  DISTRIBUTION  STATEMENT  (ot  Ihlm  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  tho  mbotrmct  ontormd  In  Block  ao,  II  dltloront  from  Report) 


It.  supplementary  NOTES 


It.  K  EY  WORDS  ( Continue  on  tovotoo  otdo  II  nocooomty  end  tdonttty  by  block  numbmr) 

nonlinear  Bernstein-Greene-Kruskal  (BGK)  wave  equilibria 
global  energy  and  momentum  conservation  constraints  ‘ 

_ •» _ 

20.  ABSTRACT  (Cantlnuo  on  «»«•••  otdo  II  NMIW|>  and  IdlUttf  br  Sloe*  itunbarj  * 

The  class  of  Bernstein-Greene-Kruskal  (BGK)  solutions  to  the 
nonlinear  Vlasov-Poisson  equations  are  examined  within  the  context  of 
the  conservation  of  (spatially  averaged)  number,  momentum  and  total 

DD  I  JAN  7}  1473  coition  OF  1  NOV  •»  is  obsolete  Unclassified 


S/N  0102-  LF-OU-SSOi 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whit  Data  M«M) 


SECURITY  CLASSIFICATION  OF  THIS  PACE  ('•Ran  Data  ent.r.tf) 


20.  Abstract,  continued 


energy,  imposed  as  ancillary  global  constraints  that  connect  the  final 
saturated  BGK  state  to  a  specified  initial  distribution  function 
f(x,v,0).  While  imposing  three  conservation  constraints  of  course  does 
not  uniquely  determine  the  final  BGK  state,  it  does  remove  a  large 
degree  of  ambiguity  as  to  whether  particular  classes  of  solutions  are 
accessible  from  given  initial  conditions.  It  also  permits  a  determina¬ 
tion  of  important  features  of  the  final  BGK  state  (e.g.,  saturation 
amplitude,  wave  phase  velocity,  etc.)  in  terms  of  properties  of  the  initial 
distribution  function.  1  The  specific  example  of  a  low-density  electron 
beam  propagating  through  a  background  plasma  with  f(x,v,0)  =  (v)  + 

nb6(v-Vb)  is  analyzed  in  detail  assuming  a  saturated  BGK  state  with 
monochromatic  waveform  $q(x')  =  $QCOSkx' ,  and  monoenergetic  distribution 

±  1/2 

of  untrapped  electrons  with  Fu(H)  =  nu(2mHQ)  6(H-Ho)  in  the  wave 

2 

frame.  Here  x'  =  x  -  (aj/k)t  and  K  =  (m/2)(v-w/k)  -  ei)>0(x')  is  the 

energy. 
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L.  INTRODUCTION  AND  SUMMARY 

The  class  of  Bernstein-Greene-Kruskal  (BGK)  solutions^  Co  the 

Vlasov-Poisson  equations  refer  to  exact  nonlinear  traveling-wave 

solutions  that  have  a  stationary  waveform  in  a  frame  of  reference  moving 

with  constant  phase  velocity  vp  =  oj/k  =  const.  The  existence  of  such 

1-4 

solutions  has  been  known  for  some  time,  and  their  properties  analyzed 
theoretically  in  various  regimes  of  interest.  Nonetheless,  nonlinear 
BGK  waves  have  remained  primarily  a  curiosity  in  the  annals  of  plasma 
physics,  particularly  because  of  the  great  generality  and  latitude  in 
specifying  the  final  BGK  state,  and  the  lack  of  connection  to  initial 
conditions  and  system  dynamics  prior  to  saturation.  The  purpose 
of  this  article  is  to  reexamine  the  class  of  BGK  solutions  to  the  nonlinear 
Vlasov-Poisson  equations  within  the  context  of  the  conservation  of 
(spatially  averaged)  number,  momentum  and  total  energy,  imposed  as 
ancillary  global  constraints  that  connect  the  final  saturated  BGK  state 
to  a  specified  initial  distribution  function  f(x,v,0).  While  imposing 
three  conservation  constraints  of  course  does  not  uniquely  determine 
the  final  BGK  state,  it  does  remove  a  large  degree  of  ambiguity  as  to 
whether  particular  classes  of  solutions  are  accessible  from  given 
initial  conditions.  It  also  permits  a  determination  of  important 
features  of  the  final  BGK  state  (e.g.,  saturation  amplitude,  wave  phase 
velocity,  etc.)  in  terms  of  properties  of  the  initial  distribution 
function. 

The  organization  of  this  paper  is  the  following.  'In  Secs.  2  and 
3,  we  review  the  general  class  of  BGK  solutions  to  the  nonlinear 
Vlasov-Poisson  equations,  including  the  conservation  of  average  number 
[Eq.  (3)],  momentum  [Eqs.  (4)  and  (25)],  and  total  energy  * [Eqs .  (5) 
and  (26)]  as  additional  global  constraint  conditions.  For  simplicity. 
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the  positive  ions  are  treated  as  an  infinitely  massive  (m.  -*■  °°) 
background  that  provides  overall  charge  neutrality.  In  Sec.  3,  for  spec¬ 
ificity,  it  is  assumed  that  the  potential  waveform  of  the  saturated 
state  is  purely  mino chroma tic  with  [Eq.  (10)] 

$0(x'  )  =  4>0coskx'  , 

where  x'  =  x-v  t,  v  =  10/k  =  const,  is  the  phase  velocity,  and  = 

p  p  Q 

const,  is  the  wave  amplitude.  In  this  regard,  it  should  be  emphasized 

that  completely  analogous  constraint  conditions  exist  for  general 

waveform  $q(x').  The  specific  choice  ^(x')  =  (fgCOskx'  is  used  to 

illustrate  the  method.  To  summarize  the  essential  results  in  Sec.  3, 

we  note  the  following.  First,  from  Eqs.  (23)  and  (24),  a  specif ication 

of  the  distribution  of  untrapped  electrons  F^(H)  =  f> (H)  +  f^ (H) 

[for  e<f>Q  <  H  <  °°]  in  the  saturated  state  uniquely  determines  the 

distribution  of  trapped  electrons  F^(H)  [for  -e4>g(x')  <  H  <  e$Q]. 

2 

Here,  H  =  (m/2)v'  -  e<j>g(x')  is  the  energy  in  the  wave  frame. 

Second,  for  the  periodic  waveform  4>g(x')  =  $gCoskx' ,  number  conservation 

follows  trivially  from  Poisson's  equation  with  =  ng  [Eq.  (21)]. 

Finally,  the  untrapped-  and  trapped-particle  distribution  functions 
+ 

Fu(H)  and  Fj(H)  are  related  to  the  initial  conditions  f(x,v,0)  through 
momentum  and  energy  conservation  [Eqs.  (23)  and  (26)].  It  is  these 
additional  constraint  equations  which  reduce  the  generality  of  the  final 
saturated  BGK  state  and  connect  properties  of  the  saturated  state  to  the 
initial  conditions. 

To  illustrate  with  a  specific  example,  in  Secs.  4  and  5  we  consider 
initial  conditions  corresponding  to  a  low-density  electron  beam 
propagating  through  a  cold  background  plasma.  For  small  initial 
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perturbations,  the  initial  electron  distribution  function  is  given  by 
[Eq.  (31)] 


f(x,v,0)  =  (nQ-nb)5(v)  +  nfe6 (v-Vfc) , 


in  the  laboratory  frame.  It  is  also  assumed  that  the  saturated  state  is 
characterized  by  monoenergetic,  backward-moving  (v*  =  v-v^  <  0)  electrons 
in  the  wave  frame  with  untrapped  distribution  function  [Eq.  (33)] 

Fu(H)  “  n0(2mH0)1/26(K-HQ), 

for  e<(>0  <  H  <  «>.  Here  n^,  n^,  V^,  nu,  and  are  constants. 

Making  use  of  Eqs.  (22)  and  (33),  the  untrapped  electron  density 
n^x'),  trapped  electron  density  n^,(x'),  and  trapped  electron 
distribution  function  FT(K)  are  exactly  determined  by  Eqs.  (34)  -  (37). 

As  a  final  assumption  (unrelated  to  the  conservation  constraints), 
we  consider  the  case  where  the  entire  beam  is  trapped  by  the 
electrostatic  potential,  i.e.,  <n-p(x')>  =  nb,  where  <...>  = 

This  relates  the  constant  in  Eq.  (33)  to  n^,  n^,  and  the  dimensionless 
parameter  e^/H^  [Eq.  (41)]. 

In  Secs.  4  and  5,  we  also  make  use  of  Eqs.  (31),  (36),  and  (41), 

and  the  expressions  for  FT(H)  [Eq.  (37)]  and  F“(H)  [Eq.  (33)],  to  simplify 

the  momentum  and  energy  constraints  in  Eqs.  (25)  and  (26).  This  leads 

to  three  equations  [Eqs.  (42),  (48),  and  (49)]  which  connect  the  four 

quantities  (saturation  amplitude),  (untrapped  electron  energy  in 

2  2 

the  wave  frame),  v  «  w/k  (wave  phase  velocity)  and  k  4>n/16 n  (average 
P  u 

electric  field  energy  density)  to  the  initial  beam  kinetic  energy 

o 

density  n^mV^/2,  beam  density  n^ ,  and  total  plasma  density  n^.  These 
three  equations  with  four  unknowns  permit  a  unique  determination  of 

2*2  a 

the  three  dimensionless  quantities  k  $q/  16TrnQHQ  =  W^,  t 


f2ll/k  d(kx') 

J0  2" 
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2  1/2 

and  vp/vo  =  (mvp/2H(p  characterizing  the  final  state  in  terms  of 
normalized  beam  density  nb/nQ  and  normalized  beam  kinetic  energy  nbmVb/2n 
^  [Eqs .  (42),  (53),  and  (54)].  For  example,  shown  in  Figs.  3  and  4, 
are  plots  of  the  Wp  =  const,  and  =  const,  universal  contours  in 

the  parameter  space  (e$0/Hfl,  nb/nQ)  obtained  numerically  from  Eqs.  (42) 
and  (53) , 


W 


F 


k2*2 
k  *0 


16™oHo 


”(1  +  e^Q/H0)1/2 
F(tt,  K)(l-e^0/H0)1/2 


1 


and 


*b  5 


vVf 

noHo 


m(l  +  ed>0/H0) 


1/2 


.  1/2 


.  _5  _0 

8  »0 


) 


F(tt  ,k)  (l-e<(iO/Ho) 

i  eVHo  _  /  _  nb^2(1+40/Ho) 

8  (l-nb/nQ)  \  nn  )  J2,_  ’ 


F  (it  ,k  ) 


2 

where  F(tt,k)  is  the  elliptic  integral  defined  in  Eq.  (40)  with  k  = 
(2ci{iq/Hq) /  (1+ciJiq/Hq)  .  If  values  of  the  dimensionless  beam  parameters 
nb^n0  and  are  Prescri^ed»  then  the  corresponding  self-consistent 

*  o 

values  of  e<p0/H0  and  Wp  follow  directly  from  Fig.  4  [Eq.  (53)]  and 

Fig.  3  [Eq.  (42)],  respectively.  For  small  beam  density  and 

potential  amplitude  with  nb/nQ,  e q / Hq  <<  1,  the  exact  expressions 

for  Wp  [Eq.  (42)]  and  [Eq.  (53)]  reduce  approximately  to  Eqs.  (57) 

and  (58),  which  can  be  used  to  estimate  (for  example)  c^q/Hq  directly 

2 

in  terms  of  nb/nQ  and  nbmVb/2ngHQ  [Eq.  (59)]. 

To  summarize,  for  the  specific  example  analyzed  in  Secs.  4  and  5, 
it  is  clear  that  the  ancillary  application  of  global  conservation 
constraints  can  provide  a  very  useful  connection  between  the  initial 
plasma  state  and  properties  of  the  final  saturated  BGK  state  (such  as 
saturation  amplitude,  average  field  energy  density,  etc.).  Of  course, 
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the  present  theory  does  not  prove  that  a  particular  saturated  BGK  state 
is  dynamically  accessible  from  assumed  initial  conditions.  However,  if 
a  particular  saturated  BGK  state  is  dynamically  accessible,  the  present 
type  of  analysis  does  remove  a  large  degree  of  the  ambiguity  associated 
with  the  final  state  and  connects  the  saturation  amplitude,  wave  phase 
velocity,  etc. ,  directly  to  initial  properties  of  the  plasma. 

We  reiterate  that  the  specific  choice  of  $q(x'),  F~(H)  and  f(x,v,0) 
were  used  to  illustrate  the  practical  utility  of  the  method,  and  other 
examples  can  be  analyzed  in  a  completely  analogous  manner. 
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2.  THEORETICAL  MODEL  AND  ASSUMPTIONS 


The  present  analysis  assumes  electrostatic  perturbations  in  a 
nonrelativistic  unmagnetized  plasma  with  one-dimensional  spatial 
variations  (3/3x  1  0).  The  ions  are  treated  as  an  infinitely 
massive  (nu  -+  <=°)  background  with  uniform  density  n^  that  provides 
overall  charge  neutrality.  In  the  two-dimensional  phase  space  (x,v) , 
the  electron  distribution  function  f(x,v,t)  evolves  according  to  the 
Vlasov  equation, 


h+  v^+  *  ML] 

^  3  t  3x  a  3x  3v  J 


f(x,v,t)  =  0 


(1) 


where  E(x,t)  =  -3$  (x,  t)/3x  is  the  electric  field ,  and  i)>(x,t)  solves 
Poisson's  equation, 


32<fr  _ 

3x2  “ 


4me  I 


dvf (x, v, t)  - 


(2) 


In  Eqs.  (1)  arid  (2),  -e  is  the  electron  charge  and  m  is  the  electron 
rest  mass. 

Equations  (1)  and  (2)  possess  three  exact  nonlinear  global 
invariants  corresponding  to  conservation  of  particle  number,  momentum, 
and  total  energy.  Assuming  that  f(x,v,t)  and  <f>(x,t)  are  characterized 
by  some  basic  spatial  periodicity  length  L,  these  conservation  relations 
can  be  expressed  as 
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dv  f  (x,v,t)  =  C1 


const .  , 


(3) 


dx 

L 


dv  v  f(x,v,t)  =  =  const. 


(A) 


'.L  dV  2  v2f  (x’v>t)  +fc(k^(x't))) 


=  const .  ,  (5) 


where  =  n„  follows  from  overall  charge  neutrality  [Eq.  (2)],  and  dx 

J  L 

extends  over  the  basic  periodicity  length  L.  However  complicated  the 
nonlinear  evolution  of  f(x,v,t)  and  4>(x,t),  the  quantities  C^,  , 

and  defined  in  Eqs.  (3)  -  (5)  remain  exactly  conserved,  and  the 
values  of  C^,  ,  and  can  be  determined  directly  in  terms  of  the 

initial  conditions  f(x,v,0)  and  (fi(x,0),  with  =  nQ  as  previously  noted. 

With  regard  to  the  nonlinear  BGK  analysis  of  Eqs.  (1)  and  (2) 
presented  in  Secs.  3-5,  we  will  find  that  the  conservation  equations 
( 3 ) — ( 5 )  provide  very  important  nonlinear  constraints  on  the  evolution 
of  the  system.  Indeed,  these  constraints  greatly  reduce  the  generality 
and  arbitrary  nature  of  the  saturated  BGK  solutions  for  specified 
initial  conditions. 
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3.  NONLINEAR  TRAVELING  WAVE  SOLUTIONS 


A.  Solutions  to  Vlasov-Poisson  Equations 

1-3 

In  the  traditional  sense  of  Bernstein,  Greene  and  Kruskal, 
we  investigate  saturated  traveling-wave  solutions  to  the  nonlinear 
Vlasov-Poisson  equations  (1)  and  (2)  that  are  stationary  in  a  frame  of 
reference  moving  with  phase  velocity 


v  =  —  =  const . 
P  k 


Defining 


x'  =  x  -  v  t  , 
P 


V '  =  v  -  V  , 
p 


t'  =  t  , 


(6) 


(7) 


and  denoting  the  stationary  solutions  (3/3t’  =  0)  in  the  moving  frame 
by  fg(x' ,v')  and  ^(x1),  it  is  readily  shown  from  Eqs .  (1)  and  (2) 
that  fgCx'.v')  and  $q(x')  satisfy 


C  3  e  3*0(X,)  3  \ 

(v'  +  =  °  ’ 


(8) 


3  *0(x') 


3x 


,2 


=  Airef  dv'  fQ(x' ,v' )  -  nJ  , 


(9) 


where  n(x')  =  i  dv'fo(x',v')  is  the  electron  density.  For  specificity, 

•  — oo 

we  restrict  the  subsequent  analysis  to  the  case  where  <$>q(x')  has 


the  waveform 


where  the  constants  and  k  are  taken  to  be  positive  without  loss 
of  generality.  That  is,  it  is  assumed  that  the  nonlinear 
saturated  state  <p^(x')  has  a  purely  monochromatic  waveform  with 
constant  amplitude  q  =  const.  It  should  be  emphasized,  however, 
that  the  formalism  and  procedures  developed  in  the  present  analysis, 
including  use  of  the  ancillary  constraint  equations  (3)  -  (5),  can  also 
be  applied  to  more  general  nonlinear  periodic  functions  $  (x'). 

The  general  nonlinear  solution  to  Eq .  (8)  can  be  expressed  as 


f0(x'v')  -  f> (H)  «(v')  +  f< (H)  «  (-v’>  , 


(11) 


where  H  is  the  energy. 


H  =  f  v-2  -  e*0(x’) 
and  «(v')  is  the  Heaviside  step  function. 


(12) 


«(v') 


+1  ,  v'  >  0  , 

0  ,  v’  <  0  . 


(13) 


For  the  electrons  trapped  in  the  electrostatic  potential 
[-e<t>g(x')  <  H  <  e<pQ ] ,  it  necessarily  follows  that  the  phase- 
space  density  of  forward-moving  (v'  >  0)  and  backward-moving  (v'<  0) 
particles  are  equal.  That  is,  f  (H)  =  f< (H)  for  the  trapped  electrons, 
and  the  phase-space  density  of  trapped  electrons  can  be  expressed  as 


fT(x’,v')  =  ~  Ft(H)[«(v’)  +  €(-v’)]  , 


(H) 


for  -e$0(x' )  <  H  <  e$Q  , 


where 
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Ft(H)  e  2f>(H)  =  2^00  ,  (15) 

for  -e<J>Q(x')  <  H  <  e$Q.  On  the  other  hand,  for  the  untrapped  electrons 
t e<J> q  <  H  <  “],  f>  (H)  and  f (H)  can  be  specified  independently,  and  the 
phase-space  density  of  untrapped  electrons  can  be  expressed  as 

fu(x',v')  =  j  F+(H)[®(v')  +  «(-v')] 

+  y  F‘(H)[«(v')  -  <B(-v'))  ,  (16) 

for  e^Q  <  H  < 


where 

F*(H)  2  f> (H)  +  f< (K)  , 
F"(H)  2  f? (H)  -  f< (H)  , 


(17) 


for  e<j> q  <  H  <  oo. 

Making  use  of  Eqs.  (14)  -  (17)  and  dv'  =  ±dH’/[2m(H  4-  e^  (x' ) )  ]  , 

it  is  readily  shown  from  Eqs.  (9)  -  (17)  that  the  Poisson  equation  (9) 
can  be  expressed  as 


-k  $ 


Q(x')  =  4iTe[nu(x' )  + 


nT(x')  - 


(18) 


where  the  untrapped-  and  trapped-electron  densities  are  defined  by 


and 


nu(x’)  = 


nT(x')  = 


dH  F+(H) 
u 


e<j>0  [2m(H  +  e<}>0(x' ))] 


1/2  ’ 


dH  Ft(H) 


-e4>Q(x')  [2m(li  +  eij>0(x' ))  ] 


1/2 


(19) 


(20) 


In  Eqs.  (18)  -  (20),  the  electrostatic  potential  is  <£q(x')  =  ^coskx'  , 
and  the  integral  constraint 


dx' 


Y“  fnu^x’)  +  nT(x')]  =  n0  =  C1 


(21) 
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follows  directly  from  Eqs.  (3),  (9),  and  (10),  where  L  =  2it/k  and 
C  r2ir/k 

dx'  =  dx'. 

•'l  0 

Note  from  Eq.  (18)  that  once  ny(x')  is  calculated  from  Eq .  (19), 
then  the  density  of  trapped  electrons  is  determined  directly  from  Eq.  (18) 
with 

k2  - 

nT(x')  -  nQ  -  nu(x')  -  ^  $Qcos  kx’ .  (22) 


Moreover,  a  particularly  important  theorem  can  be  proved  from  Eqs. 

(18)  -  (20).  For  specified  untrapped  distribution  F+(H)  = 

f >  (H)  +  f<(H)[e<}>Q  <  H  <  <=],  it  is  straightforward  to  show  that  the 

2 

distribution  of  trapped  electrons  is  uniquely  determined  from 


0  =  g(-e4>Q)  =  nT[kx'  =  (2n+l)ir] 


n0  +  4ne 


k2^0  |*  dH  F*(H) 


e$Q  [2m(H-e4>0)] 


1/2  ’ 


where  n  =  0,  1,  2,  ...  ,  is  an  integer. 


To  summarize,  for  specified  untrapped  distribution  F^ (H)  =  f> (K)  + 
f  (H)  [ebQ  <  H  <  “],  and  monochromatic  waveform  (x ’ )  =  (j^coskx' , 
the  trapped-electron  distribution  function  F  (H)  is  determined  from 


Eq.  (23),  subject  to  the  constraint  on  wavenumber  k,  amplitude 
and  untrapped  distribution  F*(H)  imposed  by  Eq.  (24). 
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B.  Constraint  Equations 


Assuming  the  saturated  BGK  equilibrium  specified  by  Eqs.  (10), 
(14),  (16),  and  (23)  is  accessible  from  given  initial  conditions 
f(x,v,0)  and  4>(x,0),  the  conservation  relations  (3)  -  (5)  provide 
important  ancillary  constraints  on  the  nonlinear  evolution  of  the 
system.  The  density  conservation  constraint  (21)  with  =  n^ 
follows  trivially  from  Eqs.  (18)  -  (20)  and  (22)  for  the  periodic 
potential  <J>q(x')  =  4>QCOskx'  .  Expressed  in  a  frame  of  reference  moving 
with  velocity  v^  =  oj/k,  the  momentum  and  energy  constraints  (4)  and  (5) 
can  be  expressed  after  some  straightforward  algebra  as 


.  dH  F  (H)  =  C2  -  Cxv  ,  (23) 

e*0 


and 


dx 

L 


,  fr40  dH  H  Ft(H) 


r«>  dH  H  F  (H) 


-e<j>0  [2m(H+e«{i0)  ] 
m  2 


1/2 


e~Q  [2m(H+e$0)]J 


,2-2 

i k  »0 

2  8tt 


(26) 


=  C,  -  C,  T  V 
3  1  2  p 


mv  (C„ 
P  2 


ciV  ■ 


where 


"0 


coskx' 


dx 

L 


2n/k 


^  ,  and  C.,  =  n„.  In  obtaining 

47T  1  u 


Eqs.  (25)  and  (26),  use  has  been  made  of  the  identities  v  =  v'  +  v  , 


mv^/2  =  mv,2/2  +  mVp/2  +  mv'v  =  H  +  mv^/2  +  mv'Vp  +  e$g(x'), 


4>q/2 ,  and 
2~  2 


Vp/i.  .  mv  Vp  -  .  “'p/  -p 

$0(x')[nT(x')  +  nu(x’)]  =  -(k2/4ne) 


dx* 
L 

dx'  2  / 

—  V*  >  - 


;(*') 


-k24-2/8Tre  [Eq.  (18)].  Moreover,  F^(H)  and  F~(H)  are  defined  in  Eqs. 
(15)  and  (17). 

To  simplify  the  subsequent  analytic  calculations,  it  is  useful  to 
define 


H'  =  H  +  e$0(x’)  ,  (27) 

and 

4>(x' )  =  d>0<xf )  +  <f>0  .  (28) 


-12- 


In  this  case,  the  conservation  of  momentum  and  energy,  Eqs.  (25)  and  (26), 
can  be  expressed  in  the  equivalent  form 


[  dx_i_  i  r 

L  m 


e$ 


dH'Fu(H’-e$0)  =  C2  -  nQvp  , 


(29) 


and 


[  dx' 

|,e<t|  dH'H'FT(H'-e$0) 

r®  dH'H'F  (H'-e^h 
u  0 

JL  L 

0  (2mH')1/2  J 

e$  (2mH')1/2 

2  8j 


(30) 


=  C3  "  n0  2  Vp  "  ffiVp(C2  "  n0Vp}' 


To  summarize,  a  specification  of  initial  conditions  f(x,v,0) 
and  <(>(x,0)  determines  the  values  of  C£  and  from  Eqs.  (4)  and  (5). 
Moreover,  for  given  untrapped  distribution  function  F“(H)  and 
monochromatic  waveform  <pQ (x’ )  =  (j>QCOskx'  ,  the  trapped-particle  distribution 
function  F^,(H)  is  uniquely  determined  from  Eqs.  (23)  and  (24).  Finally, 
the  distributions  F^(H)  and  F“(H)  in  the  saturated  state  are  further 
connected  to  the  initial  conditions  through  the  momentum  and  energy 
constraint  equations  (29)  and  (30).  As  will  become  evident  in  Sec.  4, 
these  constraint  equations  greatly  reduce  the  generality  and  arbitrariness 
of  the  final  saturated  state  for  given  initial  conditions. 
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4.  ANALYTIC  EXAMPLE  OF  SATURATED  BGK  EQUILIBRIUM 


A.  Initial  Conditions  and  Saturated  State 


It  is  clear  that  the  formalism  developed  in  Secs.  2  and  3  can  be 
applied  to  a  wide  variety  of  initial  conditions  f(x,v,0)  and  saturated 
states  fQ(x',v')  subject  to  the  constraint  equations  (29)  and  (30) 
and  the  assumption  of  monochromatic  waveform  4>q(x')  =  ^Qcoskx'  [Eq.  (10)] 
For  purposes  of  illustration,  we  consider  the  example  of  a  low-density 
(n^  <  Oq)  electron  beam  propagating  through  a  cold  background  plasma. 

For  small  initial  perturbations,  the  initial  electron  distribution 
function  is  modeled  by  [Fig.  1], 


f (x, v , t=0)  =  (n0-nb)6(v)  +  nb6(v-Vb)  , 


(31) 


in  the  laboratory  frame.  From  Eqs.  (4),  (5),  and  (31),  the  constants 
C2  and  C^  are  given  by 


C2  “ 


Vb 


C 


3 


nb  f  V 


2 

b  ’ 


(32) 


for  the  example  considered  here. 

For  present  purposes,  it  is  also  assumed  that  the  saturated 
state  is  characterized  by  monoenerget ic  backward-moving  (v'<  0) 

1/2 

electrons  in  the  wave  frame  with  f?(H)  =  0  and  f^(H)  =  nu(2mHQ)  6(H-Hq) 

2 

for  e^Q  <  H  <  ».  Here,  ny  and  Hq  =  mV^/2  are  positive  constants, 
and  from  Eq.  (17)  the  untrapped  distribution  function  F~(H)  can  be 
expressed  as 

F+(H)  =  F~ (H)  =  n  (2mHf,)1/25(H-H_),  (33) 

u  u  u  U  U 

for  edg  <  H  < 
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Making  use  of  Eqs.  (19)  and  (33),  the  density  of  untrapped  electrons  is 
given  by 


n  (x’) 
u 


[l+e$0(x')/K0J 


1/2  ’ 


(34) 


and  from  Eqs.  (22)  and  (34)  the  density  of  trapped  electrons  is  given 
by 


nT(x,)  =  n0  '  4^  V°Skx’  ‘ 


[  l-f-e<|>  0  (x  *  ) /Hq  ] 


1/2 


(35) 


where  <t>g(x')  =  <t>gCoskx' .  For  the  choice  of  untrapped  distribution 
function  F+(H)  in  Eq.  (33),  the  constraint  condition  (24)  can  be 
expressed  as 


4tre 


n 

u 


which  is  equivalent  to  nT[kx'  =  (2n+l)iT]  =  0  for  n  =  0,  ±1,  ±2, 
[Eq.  (35)]. 


(36) 


Finally,  for  the  choice  of  F~(H)  in  Eq.  (33),  it  can  be  shown 
after  some  straightforward  algebra  that  the  corresponding  self- 
consistent  distribution  of  trapped  electrons  F^(H)  is  given  by  [Eq.  (23)] 


Ft(H)  .  S&lH •  <*V“>1/2 


,2  n 

k  +  u 


27re2  (H0“H)  (l-e$0/H0)1/2 


(37) 


for  -e4>0(x' )  <  H  <  e4)Q. 


Note  from  Eq.  (37)  that  F^(H“e$Q)  *  0  [Fig.  2].  Moreover,  the  quantity 

2  - 

k  occurring  in  Eq.  (37)  can  be  eliminated  in  favor  of  $q,  n^,  ng,  etc.  , 

by  making  use  of  Eq .  (36). 
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I 


From  Eq. 

dx'  ,  , , 
—  nT(x  ) 


(35),  the  average  density  of  trapped  electrons  <nT(x')>  = 
can  be  expressed  as 


<nT(x')>  =  nQ  -  <nu(x')> 


n  _  '  —  F(lT  ,K  ) 

0  "u  7T  -  .1/2  * 

(l+e$0/H0) 


(38) 


where 

2  2e*0/H0 

k  =  - - -  ,  (39) 

<1+eW 


and  F(tt,k:)  is  the  elliptic  integral  defined  by 


F(tt,k) 


da 


„  ,.2.2  ,1/2  ‘ 

0  (1-k  sin  a) 


(40) 


If  we  impose  the  additional  constraint  that  the  entire  electron  beam  is 
trapped  by  the  electrostatic  wave,  i.e.,  <nT(x')>  =  nb ,  then  it  follows 
from  Eq.  (38)  that 


1  F(tt,<) 


U  TT 


(1+eW 


1/2 


n0 


(41) 


which  expresses  n^  directly  in  terms  of  nQ>n^  and  c^q/Hq. 

For  present  purposes,  we  assume  that  Eq.  (41)  is  satisfied, 
which  is  equivalent  to  the  requirement  that  the  entire  beam  is  trapped 
by  the  electrostatic  wave.  In  this  regard,  it  should  be  emphasized 
that  Eq.  (41)  is  an  ancillary  constraint  imposed  in  addition  to  the  con¬ 
servation  of  momentum  and  energy  in  Eqs.  (29)  and  (30).  For  future 
reference,  it  is  convenient  to  eliminate  the  constant  from 
Eqs.  (36)  and  (41).  Defining  Wp  i  k2^/  16irn0H0 ,  which  is  a  dimensionless 


| 
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measure  of  the  average  electric  field  energy  density,  it  follows  from 
Eqs.  (36)  and  (41)  that 


u2:2 

'  k  *0 
WF  '  16irn.H 


0  0 


i  f!o 

A  Hft 


7T(l+e4>0/H0) 


1/2 


F(tt  ,k  )  ( 1—  eq> q/ Hq ) 


1/2 


-  1 


(42) 


From  Eq.  (42),  it  is  straightforward  to  determine  numerically  [Fig.  3] 
the  allowed  W  =const.  contours  in  the  parameter  space  (e^/H.,  n, /n»). 

With  regard  to  Fig.  3,  the  Wp=const.  contours  are  plotted  for  parameters 
spanning  the  range  0  <  ^q/Hq  <  1  and  0  <  n^/n^  <  1-  (Keep  in  mind 
that  >  0  has  been  assumed  without  loss  of  generality.)  Moreover,  since 
Wj,  >  0  is  required  for  physically  allowed  solutions,  it  follows  from 
Eq.  (42)  that 


nb  F(7r,K)(l-e$0/H0)1/2 

°0  "  rrd+e^/Ho)172 


(43) 


In  Sec.  4.B,  we  will  find  that  energy  and  momentum  conservation  will 
further  restrict  the  allowed  values  of  (W^,  b$q/Hq,  n^/n^)  relative 
to  the  constraint  in  Eq.  (42). 


B.  Energy  and  Momentum  Conservation 

In  this  section,  we  impose  the  momentum  and  energy  conservation 

constraints  in  Eqs.  (29)  and  (30)  for  the  choice  of  self-consistent 

distributions  F^(H)  and  F,j(H)  in  Eqs.  (37)  and  (33).  Substituting 

Eqs.  (33)  and  (37)  into  Eqs.  (29)  and  (30),  and  making  use  of  C 2  =  n^V^ 
2 

and  C.  =  n,mV,/2  for  the  choice  of  initial  conditions  in  Eq.  (31), 

J  b  b 

it  follows  after  some  straightforward  algebraic  manipulation  that 
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/__0  \ 

-  {—)  nu  =  nbvb  '  n0Vp  > 


„  k2  /I 

n0H0  8tt  2  \  2 

2iTe 


(  2  H0e$0  16  6  ^0  )  2  e*0  j  ,  1/2 

(l-e$  /H  ) 

,  (45) 


,2-2 

■  1  0  m2  m  .,2 

+  2lT  =  n02Vp  +  nb2Vb'  nbmVpVb 


The  first  square-bracket  term  in  Eq.  (45)  represents  the  contribution 

from  the  untrapped  electrons  in  Eq.  (30),  whereas  the  second  square- 

bracket  term  represents  the  contribution  from  the  trapped  electrons. 

1/2 

Defining  =  (2Hg/m)  ,  and  making  use  of  Eq.  (36)  to  eliminate 

2'-  2 

k  4>q  in  the  term  proportional  to  k  H^ei^g  in  Eq.  (45),  the  conservation 
equations  (44)  and  (45)  can  be  expressed  in  the  equivalent  form 


-n  Vn  =  n,  V.  -  nnv  ,  (2 

u  0  b  b  0  p 

(1-eV  V  k  *0  m  2  m  2  .. 

n..H-  - r-TTT  -  —  =  nn  r  v  +  n,  —  V  -  n,mv  •  (‘ 

u  0  u-e»0/V  p  p 

If  we  further  make  use  of  Eq.  (41)  to  eliminate  n^  in  favor  of  n^.n^ 
and  e<f>g/HQ,  the  conservation  equations  (46)  and  (47)  become 


n (1+e^g/Hg) 


-(n0-VV0  - FOrTic) -  =  nbVb  -  n0Vp  ’ 


(no'nb)Ho 


*  '  1/2  2*2 
(l-e*0/2H0)  n(l+e*0/H0)A/^  k 


(l-e^/H0) 


1/2  F(h,k) 


m2,  m  ,,2  ,, 

=  n0  2  VP  +  nb  2  Vb  -  nbmvPVb 


To  summarize,  Eqs .  (42),  (48),  and  (49)  constitute  the  final 

equations  of  the  present  theory.  These  equations  connect  the  saturation 

2^2 

amplitude  average  electric  field  energy  density  k  c^/ 16* ,  untrapped 
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electron  energy  HQ  in  the  wave  frame,  and  wave  phase  velocity  v  =  u/k 
to  the  beam  density  nb>  total  plasma  density  n0>  and  the  initial  beam 

kinetic  energy  density  n,  mV  /2, 

b  b 
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5.  ANALYSIS  OF  CONSTRAINT  EQUATIONS 


A.  General  Analysis 


As  indicated  at  the  end  of  Sec.  4,  the  three  constraint  equations 

*  2 *  2 

(42),  (48),  and  (49)  relate  the  four  unknown  quantities  4>q,  k  $q/167t, 

H.,  and  v  associated  with  the  final  saturated  state  to  the 

0  p 

2 

quantities  n„,  n,  and  n.mV.  /2  that  characterize  the  initial  distribution 

U  D  D  D 

of  electrons.  Three  equations  with  four  unknowns  permit  a  unique 

-  2*2 

determination  of  the  three  dimensionless  quantities  c(}>q/Hq,  k  4>Q/167tnQHQ 

2  1/2 

and  v  /V.  =  (mv  /2H_)  .  Defining  the  dimensionless  parameters 

p  0  p  0 


n,  V  v 

*  _  b  *  _  b  *  p 

b  nQ  b  v0  P  v0 


where  Vq  e  (2Hq/iti)  ,  and  making  use  of  Eq.  (42)  to  eliminate 
2"  2 

the  k  $q/32tt  term  in  Eq.  (49),  the  momentum  and  energy  conservation 
equations  (48)  and  (49)  can  be  expressed  as 


(i-: 


nb\Tr(1+e^0/H0) 
TlJ  F(7T,K) 


=  VP  ‘  Vb  * 


l,_\\  /  _  5  %  \  ,1*0 

'  n0  F(»,K)(l-e$0/H0)1/2  8  H0  8  H0 


=  -  2n,  v  V  +  n,  V^. 

p  b  p  b  b  b 


Squaring  Eq.  (51)  and  subtracting  from  Eq.  (52),  we  eliminate 


v  -  2n,  v  V.  and  obtain  a  closed  expression  for  the  normalized  beam 
p  b  p  b 

*  2  - 

kinetic  energy  density  =  n^^  tenns  e^c/^0  an<^  nb//n0' 


This 


gives 


2  1/2 
tynVb/2  i.(l^0/H0)X^  _  5 

^  n0H0  F(tt  ,  k)  (l-e$0/HQ) 1;^2  8  H0 

i  eVHo  /,  ^  71 2(1+6 VV 


+  _ (1  _  'Ji.) 

+  8  (1-iWrO  V1  nJ 


“b'“0'  ‘  “0  F  (it  ,  <) 

Figure  4  illustrates  the  =  const,  contours  obtained  numerically 

from  Eq.  (53)  in  the  parameter  space  (e^/H^,  n^/ng).  The  universal 

curves  generated  in  Fig.  4  have  the  following  practical  use.  For 

specified  values  of  n^/n^  and  =  n^mV^/2nQHQ,  the  corresponding 

self-consistent  value  of  normalized  wave  amplitude  ®$q/Hq  can  be 

determined  directly  from  Fig.  4.  For  example,  for  =  0.1  and  n^/n^  = 

0.08,  then  =  0.25  follows  from  Fig.  4  and  Eq.  (53). 

Moreover,  for  consistent  values  of  nb/nQ,  etf’g/Hg  and  ,  Eq.  (42) 

can  then  be  used  to  calculate  the  corresponding  value  of  normalized 
-  2~  2 

field  energy  Wp  =  k  <j>g/16TrngHQ .  For  example,  for  (K^,  n^/nQ, 
ei0/H0)  =  (0.1,  0.08,  0.25),  then  Wp  =  k2p2/16TTn0H0  =  0.03  follows 
directly  from  Eq.  (42). 

Finally,  the  normalized  phase  velocity  Vp  =  v^/Vq  can  be  calculated 
directly  from  Eq.  (51)  by  making  use  of  Eq.  (53)  to  eliminate  V^. 

We  find  from  Eqs.  (51)  and  (53)  that 


f  2  v1/2  1/2 

*  <r)  (V1/2+(1--) 

p  \2Ho  /  v  V  13  v  V 


nn/  F (n ,<) 


where  ^(n^/ng,  e$g/Hg)  Is  defined  in  Eq.  (53).  For  consistent  values 
of  nfe/nQ  and  determined  from  Eq.  (53),  we  then  utilize 

Eq.  (54)  to  determine  the  corresponding  value  of  v^.  For  example, 
for  (K^,  nfe/nQ,  e^.  0/HQ)  *  (0.1,  0.08,  0.25  ),  it  follows  from  Eq .  (54) 


that  v 


0.998 
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B.  Small-Amplitude  Analysis 


In  Sec.  5. A,  we  considered  the  constraint  equations  (42),  (51), 
and  (52)  for  general  range  of  parameters  Rb  =  nbmV2/2n0Hg,  e$0/HQ 
and  n^/nQ.  In  this  section,  we  consider  the  case  where  the  beam 
density  and  energy  are  sufficiently  small  that  it  is  valid 
to  approximate 


n0  ’  H0 


«  1  , 


in  the  constraint  equations.  In  this  case,  correct  to  quadratic 
order  in  (e<f>g/Hg) ,  we  approximate 


If 

e$g 

\  , 

3/ 

e*o\ 

2  \ 

Ho 

)  + 

16  V 

h0  ;  ’ 

1  + 

1 

ebg 

1  /e*o 

2 

Ho 

-  8Ug 

<w„0r1/2.1  +  i(^H(^)  . 

in  Eqs.  (42),  and  (51)  -  (53).  Retaining  terms  to  order 
2  - 

(ei|>0/H0'>  an<^  nb^e^0^0^’  We  from  Eqs  ■  (42)  and  (53)  (for  example) 

that 


/!!o \  a. 

re*o  \ 

•  -l 

Ug  I  (2l 

.k0  ; 

na 

n.mV,  /2  n. 


I| 

r  e<>0  \ 

Ke_l0' 

,  n, 

|  ,  I  _b| 

4  ' 

'  Hq  > 

l'Ho  ' 

l+2no) 

If  and  nb/nQ  are  viewed  as  prescribed  quantities,  then  ef  /Hn  is 
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determined  from  Eq.  (58)  to  be 


4( 


+  64<v-v>1/2 


(59) 


where  n^  =  n^/ng,  and  >  n^  is 
The  expression  for  ecf>0/ HQ  in  Eq . 
(57)  to  det  ermine  the  normalized 


required  from  Eq.  (58)  for  c^/Hq  >  0. 

(55)  can  then  be  substituted  into  Eq. 

field  energy  W  . 

F 
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6.  CONCLUSIONS 


In  Secs.  2  and  3,  we  reviewed  the  general  class  of  BGK  solutions 
to  the  nonlinear  Vlasov-Poisson  equations,  including  the  conservation 
of  average  number  [Eq.  (3)],  momentum  [Eqs.  (4)  and  (25)],  and  total 
energy  [Eqs.  (5)  and  (26)]  as  additional  global  constraint  conditions. 

For  specificity,  it  was  assumed  that  the  potential  waveform  of  the 
saturated  state  is  purely  monochromatic  with  [Eq.  (10)] 

4>q(x'  )  =  4>0coskx'  , 

where  x'  =  x  -  v^t,  Vp  =  <*>/k  =  const,  is  the  phase  velocity,  and  ^  =  const, 
is  the  wave  amplitude.  To  summarize  the  essential  results  in  Sec.  3, 
we  note  the  following.  First,  from  Eqs.  (23)  and  (24),  a  specification 
of  the  distribution  of  untrapped  electrons  F+(K)  =  f>(H)  +  f< (H) 

(for  e<j>g  <  H  <  “>)  in  the  saturated  state  uniquely  determines  the  distribu¬ 
tion  of  trapped  electrons  F^,(H)  (for  -e$Q(x')  <  K  <  e<t>Q).  [Here,  H  = 

2 

(m/2)v'  -  e<(iQ(x')  is  the  energy  in  the  wave  frame.]  Second,  for  the 

periodic  waveform  4>q(x')  =  ^Qcoskx’ ,  number  conservation  follows 

trivially  from  Poisson's  equation  with  C^  =  n^  [Eq.  (21)].  Finally, 

+ 

the  untrapped-  and  trapped-par tide  distribution  functions  F^tH)  and 
F^,(H)  are  related  to  the  initial  conditions  f(x,v,0)  through 
momentum  and  energy  conservation  [Eqs.  (25)  and  (26)].  It  is  these 
additional  constraint  equations  which  reduce  the  generality  of  the 
final  saturated  BGK  state  and  connect  properties  of  the  saturated  state 
to  the  initial  conditions. 

To  illustrate  with  a  specific  example,  in  Secs.  4  and  5  we  considered 
initial  conditions  corresponding  to  a  low-density  electron  beam 
propagating  through  a  cold  background  plasma  [Eq.  (31)].  It  was  also 
assumed  that  the  saturated  state  is  characterized  by  monoenergetic , 
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backward-moving  electrons  in  the  wave  frame  [Eq.  (33)].  Making  use  .of 

Eqs.  (22)  and  (33),  the  untrapped  electron  density  n^x'),  trapped 

electron  density  nT(x'),  and  trapped  electron  distribution  function 

Ft(H)  were  exactly  determined  by  Eqs.  (34)  -  (37).  As  a  final  assumption 

(unrelated  to  the  conservation  constraints),  we  considered  the  case  where 

the  entire  beam  is  trapped  by  the  electrostatic  potential,  i.e.,  <n  (x')>  =  n 

r  2tt /k  ,  T  b’ 

where  <...>  =  X— ...  [Eq.  (41)]. 

Jq211 

In  Secs.  4  and  5,  we  also  made  use  of  Eqs.  (31),  (36),  and  (41), 
and  the  expressions  for  F^,(H)  [Eq.  (37)]  and  F“(H)  [Eq.  (33)],  to  simplify 
the  momentum  and  energy  constraints  in  Eqs.  (25)  and  (26).  This  led 
to  three  equations  Eqs.  (42),  (48),  and  (49)]  which  connect  the  four 


quantities  <j>0  (saturation  amplitude),  Kq  (untrapped  electron  energy  in 

2-  2 

the  wave  frame),  v^  =  ui/k  (wave  phase  velocity)  and  k  <t>g/16n  (average 

electric  field  energy  density)  to  the  initial  beam  kinetic  energy  density 
2 

vv2-  b  earn  density  n^ ,  and  total  plasma  density  n^.  These  three 

equations  with  four  unknowns  permitted  a  unique  determination  of  the 

-  2*2 

three  dimensionless  quantities  £<Pq/ Kq,  k  $Q/16iTnQHQ  i  W^,  and  v  /Vq  = 

2  1/2 

(mv  /2H„)  characterizing  the  final  state  in  terms  of  normalized 

p  U 

beam  density  nfa/n 

[Eqs.  (42),  (53),  and  (54)].  If  values  of  the  dimensionless  beam  parameters 
n^/n^  and  are  prescribed,  then  the  corresponding  self-consistent 
values  of  6<J>q/Hq  and  follow  directly  from  Fig.  4  [Eq.  (53)]  and  Fig.  3 
[Eq.  (42)],  respectively. 

For  the  specific  example  analyzed  in  Secs.  4  and  5,  it  is  clear 
that  the  ancillary  application  of  global  conservation  constraints  can 
provide  a  very  useful  connection  between  the  initial  plasma  state  and 
properties  of  the  final  saturated  BGK  state  (such  as  saturation  amplitude, 
average,  field  energy  density,  etc.).  Of  course,  the  present  theory 


q  and  normalized  beam  kinetic  energy  n^mV^/2nQHQ  i 
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does  not  prove  that  a  particular  saturated  BGK  state  is  dynamically 
accessible  from  assumed  initial  conditions.  However,  if  a  particular 
saturated  BGK  state  is  dynamically  accessible,  then  the  present  type  of 
analysis  does  remove  a  large  degree  of  the  ambiguity  associated  with 
the  final  state  and  connects  the  saturation  amplitude,  wave  phase  velo¬ 
city,  etc. ,  directly  to  initial  properties  of  the  plasma. 

+ 

We  reiterate  that  the  specific  choice  of  ( x  * ) ,  F  (H)  and  f(x,v,0) 

were  used  to  illustrate  the  practical  utility  of  the  method,  and  other 
examples  can  be  analyzed  in  a  completely  analogous  manner. 
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Fig.  1  Initial  beam  distribution  function  in  the  laboratory  frame 


[Eq.  (31)]. 
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2  Trapped  electron  distribution  function  in  the  wave  frame 


[Eq.  (37)]. 
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Plot  of  Wj,  «■  const,  universal  contours  in  the  parameter  space 
(e^0/H0,  nb/nQ)  [Eq.  (42)]. 


